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Abstract — This paper investigates the nature of the development 
of two-dimensional laminar flow of an incompressible fluid at the 
reversed stagnation-point. ". In this study, we revisit the problem 
of reversed stagnation-point flow over a flat plate. Proudman and 
Johnson (1962) first studied the flow and obtained an asymptotic 
solution by neglecting the viscous terms. This is no true in neglecting 
the viscous terms within the total flow field. In particular it is pointed 
out that for a plate impulsively accelerated from rest to a constant 
velocity VO that a similarity solution to the self- similar ODE is 
obtained which is noteworthy completely analytical. 

Keywords — ^reversed stagnation-point flow, similarity solutions, 
analytical solution, numerical solution 



I. Introduction 

THE full Navier-Stokes equations are difficult or impos- 
sible to obtain an exact solution in almost every real 
situation because of the analytic difficulties associated with the 
nonlinearity due to convective acceleration. The existence of 
exact solutions are fundamental not only in their own right as 
solutions of particular flows, but also are agreeable in accuracy 
checks for numerical solutions. 
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Fig. 1. Oil recovery industry 

In some simplified cases, such as a fluid travels through a 
rigid body (e.g., missile, sports ball, automobile, spaceflight 
vehicle), or in oil recovery industry crude oil that can be 
extracted from an oil field is achieved by gas injection, as 
shown in Eig. ([T]), or equivalently, an external flow impinges 
on a stationary point called stagnation-point that is on the 
surface of a submerged body in a flow, of which the velocity 
at the surface of the submerged object is zero. A stagnation- 
point flow develops and the streamline is perpendicular to 
the surface of the rigid body. The flow in the vicinity of 
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this stagnation point is characterized by Navier-Stokes equa- 
tions. By introducing coordinate variable transformation, the 
number of independent variables is reduced by one or more. 
The governing equations can be simplified to the non-linear 
ordinary differential equations and are analytic solvable. The 
classic problems of two-dimensional stagnation-point flows 
can be analyzed exactly by Hiemenz fT], one of Prandtl's first 
students. The result is an exact solution for flow directed per- 
pendicular to an infinite flat plate. Howarth 121 and Davey ISll 
extended the two-dimensional and axisymmetric flows to three 
dimensions, which is based on boundary layer approximation 
in the direction normal to the plane. 
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Fig. 2. Sketch of forward stagnation point and backward stagnation point 

If the rear of the rigid body is not tapered, a stagnation-point 
flow also develops in the rear of the body, as shown in Fig. (|2]). 
The flow in the vicinity of this reversed stagnation point is 
governed by boundary-layer separation and vortex generation 
and the reverse stagnation point flow develops. The main 
difference between these two flows is the change in the flow 
direction. Reversed stagnation-point flows against an infinite 
flat wall do not have analytic solution in two dimensions, but 
certain reverse flows have solution in three dimensions O. 

Proudman and Johnson first suggested that the convec- 
tion terms dominate in considering the inviscid equation in the 
body of the fluid. By introducing a very simple function of a 
particular similarity variable and neglecting the viscous forces 
away from the plane, they obtained an asymptotic solution 
in reversed stagnation-point flow, describing the development 
of the region of separated flow for large time t. In their 
solution, the phenomenon of separation is described near a 
plane that represents the rear- stagnation point of a cylinder 
is set in motion impulsively with a constant velocity normal 
to the surface of the plane. Robins and Howarth |5 | have 
recently extended the asymptotic solution, finding the higher 
order terms by singular perturbation methods. They indicated 
that the viscous forces cannot be ignored in the governing 
equation because of a consistent asymptotic expansion in both 



this outer inviscid region and also in the inner region near the 
plane. Smith f6l generalized the solution of Proudman and 
Johnson with both viscous and convection terms in balance 
by considering the monotonic potential flow when the time is 
relatively large. 

Numerical simulation of reversed stagnation-point flow with 
full Navier-Stokes equations has been studied in |7|. In the 
present study, the unsteady reversed stagnation-point flow is 
investigated. The flow is started impulsively in motion with 
a constant velocity away from near the stagnation point. A 
similarity solution of full Navier-Stokes equations is solved 
by applying numerical method. 

II. Flow Analysis Model 

The viscous fluid flows in a rectangular Cartesian coordi- 
nates (x, z). Fig. |3j which illustrates the motion of external 
flow directly moves perpendicular out of an infinite flat plane 
wall. The origin is the so-called stagnation point and z is the 
normal to the plane. 
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Fig. 3. Coordinate system of reversed stagnation-point flow 

By conservation of mass principle with constant physical 
properties , the equation of continuity is: 



du dv ^ 
dx dy 



(1) 



We consider the two-dimensional reversed stagnation-point 
flow in unsteady state and the flow is bounded by an infinite 
plane = 0, the fluid remains at rest when time t < 0. At 
t = 0, it starts impulsively in motion which is determined by 
the stream function 

tjj = —axy (2) 

At large distances far above the planar boundary, the ex- 
istence of the potential flow implies an inviscid boundary 
condition. It is given by 
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where u and v are the components of flow velocity, A is a 
constant proportional to Vq/L, Vq is the external flow velocity 
removing from the plane and L is the characteristic length. 



We have u = at x = and v = at y = 0, but the no-slip 
boundary at wall {y = 0) cannot be satisfied. 

For a viscid fluid the stream function, since the flow motion 
is determined by only two factors, the kinematic viscosity u 
and a, we consider the following modified stream function 
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(4a) 

(4b) 
(4c) 

where rj is the non-dimensional distance from wall and r 
is the non-dimensional time. Noting that the stream function 
automatically satisfies equation of continuity ([T]) . The Navier- 
Stokes equations 1 8 1 governing the unsteady flow with constant 
physical properties are 
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where u and v are the velocity components along x and y 
axes, and p is the density. 

Substituting u and v into the governing equations results a 
simplified partial differential equation. From the definition of 
stream function, we have 

d^l) 
dy 

d^l) 

dx 

The governing equations can be simplified by a similarity 
transformation when several independent variables appear in 
specific combinations, in flow geometries involving infinite or 
semi-infinite surfaces. This leads to rescaling, or the introduc- 
tion of dimensionless variables, converting the original sys- 
tems of partial differential equations into a partial differential 
equation. 

pdx 

(7a) 



pdy 

The pressure gradient can be again reduced by a further 



differentiation Eq. (7b) with respect to x. That is 

dxdy 

and Eq. ( [7a| reduces to 

[fr]T ~ iff]) + f firn ~ frjrjrj\rj — 0- 

The initial and boundary conditions are 

/(r,,0)=r; (r,7^0) 
/(0,t)=/,(0,t)=0 {t^O) 
/(oo,t) ~ 



(8) 



(9) 
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The last condition reduces the above differential Eq. ([9]) to 
the form 

fr]T ~ ifr]) ~^ f fr]r] ~ f-q-qr] — (H) 

with the boundary conditions 

/(0,r) = ^(0,r)=0 (12a) 
^(oo,r) = l. (12b) 

Eq. ( pTj ) is the similarity equation of the full Navier- Stokes 
equations at two-dimension reversed stagnation point. The 
coordinates x and y are vanished, leaving only a dimensionless 
variable r]. Under the boundary conditions /^(oo,r) = 1, 
when the flow is in steady state such that frjr = 0, the 
differential equation has no solution. 

III. Similarity Analysis 
A. Asymptotic solution 

When r is relatively small, Proudman and Johnson IH first 
considered the early stages of the diffusion of the initial vortex 
sheet at = 0. They suggested that, when the flow is near the 
plane region, the viscous forces are dominant, and the viscous 
term in the governing Navier-Stokes equations is important 
only near the boundary. 

On the contrary, the viscous forces were neglected away 
from the wall. The convection terms dominate the motion of 
external flow in considering the inviscid equation in the fluid. 
They considered the similarity of the inviscid equation 



fvr - ifvf + ffvv + 1 = 0. 



(13) 



Proudman and Johnson obtained a similarity solution of (13) 
is in the form 

/(r,,T)=e^F(7) (14) 
J^(7) =7- (15) 



where 7 = 



V 



Mr) 



and c is a constant of integration. Robins 



and Howarth [5J estimated the value of c to be 3.51. This 
solution describes the flow in the outer region, moving away 
from the plane with a constant velocity. It can be checked that 
the viscous term /^^^ is still small compared to the convective 
terms, so that their assumsion of neglecting the viscous term 
is still valid. 

In the inner region, the viscous term cannot be neglected 
and the no-slip condition must be satisfied on the wall. When 



r ^ 00 and r]/e^ is relatively small, the solution (15) yields 

F 



-7 



-r]e 



and 



/' = -i 



(16) 



f = -v, 

Substituting in ( pT) becomes 

f"-ff" + {f?-i 

/(O) = /'(o) 

/'(^) = -1 

This is exactly the classic stagnation-point problem 
(Hiemenz (H) by changing the sign in /. It is a third- 
order nonlinear ordinary differential equation and does not 



(17) 



have an analytic solution, and thus it is necessary to solve it 
numerically. The general features of the predicted streamline 
are sketched in Fig. (|4]). 
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Fig. 4. Streamlines of reversed stagnation-point flow 

Although an asymptotic solution was obtained, it can easily 
been observed that this is no true in neglecting the viscous 
terms within the total flow field. No exact solutions in both 
outer and inner regions were discovered. 

B. Particular Solution 

In our two-dimensional model, the fluid remains at rest 
when time t < and is set in motion at t > such that 
at large distances far above the planar boundary the potential 
flow is a constant Vq for all value of t. Both Proudman and 
Johnson |4 |, and Robins and Howarth (5] have set Vb = 1 and 
the corresponding boundary condition /^(oo,r) = 1. 

When the flow is in steady state such that frjr = 0, it was 
proven that the similarity velocity fr]{r]) cannot ultimately 
approach to 1. The differential equation has no solution. In 
this chapter, if the potential flow Vb is restricted not to be a 
constant, the boundary condition /^((X),r) results in a time 
dependent function and then we obtain another approach of 
similarity solution in reversed stagnation-point flow. 

As with the governing equation of reversed stagnation-point 
flow, we can write the stream function as 
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where A is a constant proportional to Vo{r)/L, Vo{r) is the 
external flow velocity removing from the plane and L is the 
characteristic length. These result in the gonvering equation 



[fr]T ifrj) ~^ f frjr] fr] 



0. 



After integration with respect to rj, we have 

friT ~ ifri) + 7/7777 ~ fr]r]r] — ~C(t)i 



(19) 



(20) 



Under the boundary conditions /^(oo,t) = 1, the value of 
C{t) should be a constant and equal to 1. If the boundary 



condition /^(oo, r) is restricted not to be a constant, t, follow- 
ing the procedure of |9|, a particular time-dependence function 
C{r) may be expressed in the form 



C(r) 



c 



(21) 



where c is an arbitrary constant. The partial differential equa- 
tion can be simplified by a similarity transformation when a 
new similarity variable is introduced. This converts the original 
partial differential equation into an ordinary differential equa- 
tion. For a time dependent function, we introduce the diffusion 
variable transformation 1.6.1 



f{v.r) = -^F{,) 



(22a) 
(22b) 



Here is the time combined nondimensional variable and 
F{(^) are the nondimensional velocity functions. Substitution 
of the similarity transformation yields an ordinary differential 
equation 



1 



^F^^ -F' - F''^ + FF" - F 



(23) 



Equation ( 23 ) is a third-order nonlinear ordinary differential 
equation and a key step in obtaining an analytical solution 
is to rearrange the equation in an autonomous differential 
equation. In mathematics, an autonomous differential equation 
is a system of ordinary differential equations which does not 
explicitly depend on the independent variable. 

In order to omitting the variable ^ in the differential equa- 
tion, it is recognized a change of variable 



-4= 



(24) 



and the equation becomes to an autonomous differential equa- 
tion 



QQ" - 2Q' - Q'^ - 



(25) 



In our analysis, P = is the dependent variable and Q is 



the independent variable. Equation (25) is rearranged as 

QP' -2P - - = -c^- 



(26) 



and the chain rule reduces equation ( [26| ) to a second-order 
ordinary differential equation 



^^dP ^ ^ d f dP 



-c+- (27) 



Equation ( [27] ) is analytically solvable that the solution might 
be expressed as a low order polynomial. It is suggested that 



P = a^bQ^dQ^ 



(28) 



and substituting into equation (26) and comparing the coeffi- 
cients in the powers of Q results in a system of linear algebraic 



equation 



4 

Mhd + 6^ + 26 + a6 = 

I2hd^ -bd = 
6d^ - = 

Solving the related algebraic equation, we have 



(29a) 

(29b) 
(29c) 
(29d) 
(29e) 
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(30) 



3 . ^ 3,1 

--, 6 = 0, c = -, d = - 
2' ' 4' 6 

Substituting the constant into equation ([28]) yields the first 

order differential equation 

3 1 
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(31) 



Equation ( [3T] ) is Riccati equation , which is any ordinary 
differential equation that is quadratic in the unknown function. 
The standard form of Riccati equation is 

Q' = RQ^^SQ^T (32) 

The solution of Riccati equation can be obtained by a change 
of dependent variable, where the dependent variable y is 
changed to q by 1 10] 

Q = -'-^ 

^ q R 

By identifying i?=^,6' = 0andT = 
variables in equation ( [3T] ) becomes 

q' Qq' 

(1 
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(33) 

\, the change of 



(34) 



so the equation (31) becomes an second-order linear differen- 
tial equation 

q" -\=0 (35) 



The general solution to this equation is 



A cosh- 
2 



q = A. costi- + stnti- (36) 
where A and B are arbitrary constants. Applying this solution 



B sinh- 
2 



in equation ( [24] ) leads to the general solution of equation ( [23 ) 

m = i;- ' . \::, c : (37) 



3A sinh^ 
A cosh^ 



- 3B cosh^ 
B sinh^ 



Application of the no- slip condition F{0) = leads to the 
determination of the constant B = 0, so the exact solution 
becomes 



- — 3 tanh- 
2 2 

Collecting results, the velocity functions become 



(38) 



(39) 



where = 



y is the non-dimensional distance from the 



plate. In view of (39), the flow far form the boundary (x, ^ 
00) becomes 

/(.,r)^-^(|-3) (40) 



We obtain a particular solution of the unsteady reversed 
stagnation-point flow. The above solution is obtained in the 
similarity framework for unsteady viscous flows. The ap- 



pearance of this positive factor in the first terms of (40) 
shows that this remote flow is directed toward the axis of 
symmetry and away from the plate. The second term in 
( [4Q| describes a uniform current directed toward the plate. 
An adverse pressure gradient near the wall region leads to 
boundary-layer separation and associated flow reversal. 

The particular solution is noteworthy in that it is completely 
analytical, but it is limited to the region far away from the wall 
in the presence of non-zero term F'{0) = —1. The no-slip 
boundary condition is not satisfied completely near the wall 
region. 

C. Numerical Solution 

Since the analytical solution does not satisfy the no-slip 
condition = 0, it is convenient to solve the similarity 
equation numerically. The similarity equation and the relevant 
boundary conditions are 



F' - + FF" - F'" = -c 
F(0) = F'(0) = 



(41) 



where c = | in order to satisfy the unsteady viscous flows in 
the outer region. 

This equation is a third-order nonlinear ordinary differential 
equation. It is convenient when solving an ODE system 
numerically to describe the problem in terms of a system of 
first-order equations in M ATL AB ifTTI . 

For example when solving an n^^ -order problem numeri- 
cally is common practice to reduce the equation to a system 
of n first-order equations. Then, by defining yi = F, ^2 = 
F' ^ ys = F" , the ODE reduces to the form 



dx 



2/2 

^3 



(42) 



The first task is to reduce the equation above to a system 
of first-order equations and define in MATLAB a function to 
return these. The relevant MATLAB expression for Eq. ( [42| ) 
would be: 

Listing 1. System of first-order equations 

function dy = s tagnation ( t , y ) 
c=3/4; 

dy = zeros (3,1); 
dy(l) = y(2) ; 
dy(2) = y(3) ; 
dy(3) = c-l/2*t*y(3)- 

y(3); 



-y(2)-y(2)*y(2)+y(l)* 



Later, we need to change the boundary value into initial 
value, because odeAb, a ode solver in MATLAB, can only 
solve the initial value problem. From Eq. ( [4T] ), we need 
to gauss the value of F"((S) such that F'{oo) = |. The 
commands written in MATLAB would be 



Listing 2. ODE solver 

function main 
x = 20; 

[T,Y] = ode45( ( 
-1.54306]) ; 
plot(T,Y(: ,1) , ' 

'-.') 



) stagnation , [0 x 
-' ,T,Y(: ,2) , ' — ' 



,[0 
,T,Y(: ,3) 



The numerical solution for two-dimensional stagnation- 
point flow is shown in Fig. ([5]). 




Fig. 5. Numerical solutions of reversed stagnation-point flow 



This solution is a similarity solution of the reversed 
stagnation-point flow over a flat plate. It describes an unsteady 
viscous flow in both outer and inner regions. A single dividing 
streamline plane separates streamlines approaching the plate 
from external flow streamlines. Though the viscous term F'" 
is still small compared to the convective terms in the outer 
region, this is no true in neglecting the viscous terms within 
the total flow field. The similarity velocity field is show in 
Fig.([6]). 



IV. Conclusion 

The foregoing study constitutes a similarity solution of 
the unsteady Navier-Stokes equations for reversed stagnation- 
point flow in the idealized case of an infinite plane boundary. 
In order to analyze the flow for non-zero values of x, it 
is required to convert the full Navier-Stokes equations. This 
problem is now being studied by applying numerical method. 
The solution is obtained in the classical similarity framework 
for unsteady viscous flows. 

In the case of numerical methods, a brief analysis of the 
new solution is discussed. When the flow was near the plane 
wall region, the viscous forces were dominant, and the viscous 
term in the governing Navier-Stokes equations was important 
only near the boundary. On the contrary, the viscous forces 
were negligible when they were away from the wall. 
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Fig. 6. Similarity velocity field as a function of <^ 
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